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PREFACE 


The theory of radiative transfer is important in 
| meteorology, for example in the interpretation of radiation 
measurements from satellites. Some of the techniques can 
also be applied to studies of neutron transport. 

The present study relates the theory of radiative 


, Seaneter to basic metheds in the theory of probability. It 


should be of interest to specialists in the theory of 
| radiative transfer and to those concerned with applications 


of probability theory to physics. 


the besis of the stochastic model of multipie 
scattering of photons, we consider the diffuse reflection and 
transmission of a parallel beam of radiation by a finite, 
plane-parallel, non-emitting and homegeneous atmosphere with 


conservative and isotropic scattering 


stochastic precess under consideration represents a homo— 
geneous stationary evolution in a Markovian manner with 


pespeet to the optical depth. | 
Pirst we derive the forward and the backward integro- | 


differential equations for the emission probability 
starting with the Laplace transform of these equations, we 


‘distributions from the Chapmaz 


obtain the S— and T-functions of 8. 


monedirectional illumination of the upper and the lower 
boundaries, depending on the optical depths +t, and 1, 

(O¢ TT) < %,)- The results obtained with the aid of the | 
forward equations reduce to those derived from the backward 
equations, because of the homogeneous optical properties of | 
the medium. Some new functional equations for the source 
‘funetions of the auxiliary equations are given. 
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ABILISTIC METHOD FOR PROBLEMS OF RADIATIVE PTRANSHER: 
MARK PERTY OF RADIATIVE TRANSFER 
AND OF NEUTRON DIFFUSION 


It is well known that the Markov property has been 

considered in detail in the field of Brownian motion [9], 
turbulent diffusion [2], and others [6]. 

| In a preceding paper [16], the author showed that : 
‘“Milne's problem in a semi-infinite, plane-parallel atmosphere | 
with isotrepic scattering can be solved, with the aid of the | 
Markovian stechastic model of multiple scattering of photons, ) 
‘by using Chand khar's idea [11], which is based on the 

: principle of invariance arising from the 
at infinity. 


| 


otic solution 


With the aid of the physical method based on the | 
‘principle of invariance, the problem of the diffuse reflection 


[1], Busbridge [7], and Chandraselch: 
‘mathematical discussion of the same problem in an inh 
: geneous medium has been given by Bellman and Kelaba [3], [4], ! 
[5], Sobolev [14], Buebridge [8], and Ueno [20]. | 
| In the present pap 


the | 


er we consider a stochastic process . 


ect | 


i 


ng the evolution in a tovian 


ia with 


and neutron diffusion. As a typical » we deal with the 


ion of a parallel beam of 


‘diffuse reflection and transmt 


‘radiation by a plane-parallel and non-emitting atmosphere of — 
‘finite optical thickness 1, — ty (0 < %) < T,) with | 
conservative and isotropic seattering, because the solutions — 
of all other transfer problems in a similar flat layer can : 
apparently be reduced to this one. In such a stochastic 
field the precess is considered as denoting the evolution of 
“the probability distribution of emission p(u; To»t.T,) with 


increasing optical depth « (0 ¢ 1) <¢1<¢1,)- Here | 
plus Torts, au is the probability of finding a continuous | 
| stochastic parameter u(t) between 4 and du at the | 
| level + where cos “y otes the inclination to the 

outward normal of the surface + = 1). It 


that P(H5 TT, 7) is a truncated function of t, i.@., | 
P(E T oT, »t,) is finite for K%StTLYH and | 
pP(styst,t,) = 0 for t<t) and t, < t- | 

In the present paper, we construct a Harter ian stochastic 
model of multiple er of photons and apsume the 


haracter of the evolution with 


respect to Tt. 
integre-~differential equation 
distributions from the Chapman-Kolmogo! 


reff equations approp- 
zg with the Laplace 


i 


‘f-Feller equations [12], we 
obtain the scattering and the transmission functions which — | 


vary with the optical 
the integral 


and tT): A complete set of 
functions derived from 
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the forward integro-dif: 
for the 3—- and T-functions derived from the backward 
because of the optical homegeneity of the medium. These 
equations shew a Markovian property of miltiple scattering of — 
classical particles, i.e., photon and neutron, under the | 
given boundary condition: 

It may be of interest to nete that eur above procedure 
is somewhat similar to that used by Feller [13] on boundaries | 
and lateral conditions for the Kolmog: 


vential equations are equal to these 


equation 
what happens in the final infinitesimal interval, in form— 

lating the backward equation we seek to show what happens in 

the initial infinitesimal interval. In other words, while 
‘the former involves differentiation with respect to 1,, the 
‘latter contains differentiation with respect to fp. 
Mathematically speaking, the forward and the backward | 
equations together can be derived from the auxiliary equations 
nination of the upper and the lower 


for monedirectional 112 


boundaries, because the emission probability distributions 
depend on the optical depths and t,- In what follows, 
| for convenience the various 
probability distribution, the auxiliary equation, the S~ and 

ur tions, and the % and Y—funections monodirect ional 
“Aliwueination of the upper and the lower aries, will be 
inthe sania 


tities, 1.e., the emission 


denoted respectively as the 


While the above Markevian stochastic model was first put 
forward by Ueno [16], the probabilistic idea in radiative | 
transfer can be traced out in the methed of invariance due to 
_Ambarsumian [1] and Chandrasekhar [10]. , ni 
of Chandrasekhar can be interpreted as the intensities of the 
thermore the %- and ¥-functions — 


prebability currents, and 
‘are equal to the probability distributions of emission from 
the bounding planes 1 = To and tT = Tp respectively. 


In later papers we shall reconside 
(of rediative transfer, allowing for the AO gi 
‘Vene {21]), nencoherent scatteria b» and other geometries than | 


Ea: a Ae Whe ad Bees os 


Asouming that the diffuse reflection and transmission of | 


} 


parallel rays by a finite flat layer with conservative and 
covian stochastic proces s 


isotropic scattering represen ai 
with respect to the optical 


: depth +, we can physically define the emission probability 


i 
i 


distribution as follows: Let p(H3t>,T,T) be the 


‘probability that a photon absorbed at level + (or 1, + tT) — 7) 
e-~emitted in the direction +4 (or — #) be 
(0 <#¢ 1) in the radiation emerging from the 


ve, let p(M3%,,%,T) be the | 


will be 


(‘Tet (or t= 1). thermo 
emergence probability of a photon absorbed at the level t 


(or 1, +%)— 1) in the direction~u (or +4) (0 <u ¢2)_ 


RM-2936 
ae 
omoge: cous case [20], we shall eall 
p(s ty, 4.7,) and p(t,» Fp) the and the upward 
emission prebability distributions, respectively. 
In a preceding paper [17], assuming the Mi , property 
of the diffuse reflection and transmission of radiation for 


to that used in the inh 


isetrepiec and coherent seattering in the conservative case, 


we derived the Kolmogoreff—Feller equation frem the Chapman- 
‘Kolmegereff equation, allewing for the leakage of the 

unding plane t= 1, due. 
eted towards the level 


probability current throug? 
to the transmission of radiation dire 
a 2.) T° 

In a manner similar to that used in probability theory 


(12], we find the Chapman-Kolmogeroff equ 


tion appropriate to. 


the present case to be 

: om | 

(2-1) P(t pst», ) -/ B(H' stp, t-At, t, P(A! At, t, aul, 

0 i 

for all values of A between 4, and +,- In Bq. (2.1) 
P(; Tt, T-A7,7,) and p(u/i';At,1,) are given by 


(2.2) Blu to, at, t,) = plasty, tat, 4, ) 


— At i pth" 3} Ty t, T))Ro (u/n" jan" , 
te) 


~a | | 
(2.3) p(u/a' sAt,t,) = Ry (m/e! ar + sewn — At 7 R(A"/a! a") 
0 | 


6 


In the above, the probability distribution function 


p(u,/u, itt, +7 ai, represents the probability that 4(t,) 
lies between 4, and K+ du, provided that u(t, ) = ie 
(O< TF C™<rS t,)- In (2.3), & 1s the Dirac delta 


function, and Ry Ros and R-functions are provided by 


X(u,t,), Ro(W/a') = gtr Yu,,), 


, (2.4) Ry (u/u') = aie 
(2.5) R(w/u') = Ry (u/a') + Ro (u/s'), 

were 

| x,t) = P(M5 toe T,) = X(H,t,) = P(Hs ty» T)>T), 
(2.6) 


¥(e,,) = p(H3 5.7.7) = Y(u,t)) om P(WST,»T eT )- 


The quantity R(H/1')dp dv represents the conditional 
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transition probability thet, given #', one finda in the | 
range (i, 4 + du) through a parametric interval dt of the. 


optical depth. 
The normalization condition for the conditional proba— 
bility R(u/u') in the conservative case is provided by 


| 1 

(2.7) uf R(w' Anja! = 1. 

: 0 

In the nonconservative case the above integral is less than 


unity. While an explicit appeal to the K-integral in the 
olive the ambiguity 


conservative case is needed to res 


the normalization cendition (2.7) is also useful for the 
resolution of the arbitrariness of the solution (see Ueno 


Letting At +0, we obtain the stochastic integro— | 


‘differential equation for p(s t,t. T,) 


Sp(u;t 


ne toy) | 
a a 2 p(H3t),*,7,) | 


(2.8) 

i at 

+ + x(u,t,) f p(H'5 +), 7,7, } oH 
6 


7. 


1 
2s + v(usz,) f P(H'5 4,7, To) “at 
2] 


Similarly, we can obtain the other ~differential 


equation for the upward emission probability distribution 
p(ust,,%, T») . We write 


1 3 
(2 9) P(H3745T,T,) = f B(m' 541. T4AT, ty )p(M/n' ;At,T))au', 
3] | 


(2-10) Blust,, T4841) = plas t,, TAT, T)) | 
: | 


1 
= ax [ p(n" st, t4ht, +, )R,(4/2" an", 


0 | 


(2.22) p(/u'sat,t9) = plu/u'sar, ty), 


in the parameter + are identical in the homegeneous case. 


As At +0, we have 


—+ X(u,%,) iis p(#'5,,%,T)) ree 


~1 Pe 
+ 5 ¥(u,«,) f p(h' 34,4, 7, ) re 
6 


riate to this case takes the form 


Next we derive the fo 
_Kelmegoreff equation 


| | | 
(2.13) p(#st,,.T,T,) -/ B(m' 5 t95%, At )p(u/n' pt ,At du’, | 
: A 


In (2.12) we put 


(2.24) Bltstys t.ty-At) = p(s ty, t,7,-At) 


,j 
+ AT f p(n" 5%,-4T, 7, T>)Ro (m/a" jan" , 
0 


(2-15) p(u/utzt,at) = O(m — Bt). 


Then, letting At 0, we get 


* # ¥(u,7,) a P(H'S Tt ,t.T 5) 


Similarly, we may start with anether type of Chapman— 
reff equation for the upward emission probability 
distribution, namely 


(2.17) P(H5 1,7, T,) -[ D(w' 5 t,-8t,7,15)p (MAb st at )du', 


(2.28) P (w/t 5,A7) = p(u/u' sat, t,)- 


Then as At +0, (2.17) becomes 


+ 3 xu.e,) f° p(t'3%,,1,%) ag : 
0 


On combining Eqs. (2.8) and (2.12) with Eqs. (2.16) and 
(2.19), respectively, we obtain. 


= — § pluitgstoty) + F X(H,4,) 


r} a 
. 7] P(M'5t5,t,7,) ae j 
fs) 
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dp(ust, t,t ) 


(2.21) 


1 af 
= $ run) [ p(M 557 7,) ad : 


Now we derive the backward 
'the forward equations (2.16) and (2.19). First we write 


equations corresponding to 


: 1 
(2.22) P(M5t5.T)T,) = 7 B(w' ;totAt, tT, 7, p(uAL' 54,7 )au', | 
) 


where p(i/u';t,4+) is equal to p(u/u' 47,7, ) in (2.3). 
Then, as At + 0, 


(2.23) on Zt p(u; To? *;) 


ni is 
-% X(u,»,) ; p(u'st, t,t) Ta : 
0 


Similarly, if we start with the Chapman-Kolmogoroff 


equation ylelded by 


i 
(2.24) p(s ,,T,T)) -/ Bust, t,t tht )p(u/n' +, At au’, 
6 


(2.25) B(u5T,,t, Todt) = p(w tt), Tt, ttt) 


+ At i p(e" st otht, 1,7) Ry (m/f) “FH 


(2.26) pluftstas) = O(p—nt), | 


then, in the limit as At +0 we get 


.1 
-$ ¥(u,t,) i p(H'3%5,7,T,) 
0 


| The combination of Eqs. (2.23) and (2.27) with Eqs. 
(2.8) and (2.12) prevides the following equations: 


= + P(M5t,>T,% >) 


i 
-$ x(u,e,) [ P(M'5t),7,%—) Ser 


+ ———+3= | 
| 


While Eq. (2.23) is equal to Eq. (3.23) given by 


‘Busbridge [8], so far as we know, Bes. (2.28) — (2.31) are 
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new. Furthermore, it is known that the above emission 
“probability distribution fulfills the following downward and | 
the upward auxiliary equations (see Ueno [10]): | 


(2.32) 2 — BY, (p(usty,t,4,)) = tee 


(2.33) {1 — XK), (p(ust,,t,4))) = Nie 2 , | 


where 1 is the identity operator and the truncated Hopf 
operator XK is 
' . 


(2.38) X,(2(t)) = {> £(t)e,(|t — +|)at. 
2 Be. 
Le) 


In (2.33) 3B, is the first exponential integral | 


(2.35) By (+) Nia o/h ae 
0 


In the conservative case, Eqs. (2.32) and (2.33) show that 
P(H3t.T,4,) eoincides with P(B5T) st ott \—T,T))- 


= st ~(tt9)/u | 


, | 
sg 


a eo ee 

t 

nt —(4,—*)/t 
(3.3) T(%5s7 sso) = ‘ P(Hy3 Tot, 7, Je at, 
| nm, —(tt,) /t 
(3.4) T(,sTosHoty) -{ 3 P(o3 TTT) Je . 
t. 
) 


dt. 
However, allowing for p(ust4>%,T) = P(Hst5, tot, 7,7, ) | 


| in the conservative case, we have 
(3-5) S(¥,4, so) = 8(4,, 495K .40), 
(3.6) Tt sTy sy) = T(T) Fo 5H9)- 


If the optical properties of the medium vary with the : 
optical depth (see Busbridge [8], Ueno [20]), Eqs. (3.5) and | 
(3.6) do not hold. 


From (2.32) and (2.33) we get | 


(3-7) X(u,t,) = 1+ $ : : 8(t).4 5,6") ABr 
* 
i) 


(3.8) ¥(u,t,) = o rte + bf * T(T),T, 34,4") SH 


és 


(3-9) S(T or TysHrty) = S(ty,T ity H)> 
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(3.10) P(t5.T, soy) = My, Tso). 


(45) /ty 
mar ai 5 laa the forward equation (2.19) by e 
(t,t) /a 
@ a 6 respectively, and integrating with respect 


to +t over (t),,), from Eqs. (3.1) — (3.10), we have 


pt 
(3.12) S(T, ,ToshsHy) = if X(H,€)X(Ho,t) 

. 

° 


- exp {(— (*, - t)(2 + i.) iat, 


(3-12) Mt), T,s4H_) = : xiu.t)tugetle at, 


J 


To 


| where 


(3-13) -X(u,t) = 2+ f° ef X(u,t)x(',t) 
oO Tt 
- exp (~— (* — te + Pr) at, 


a, f° X(n,t)¥(u',t)e ‘ Ms dt 


¢ 


(3.14) 


—( tt) /4 
Multiply the forward equation (2.16) by e ® 


and e = oe grate with respect to | 


= over (Tos7,)- 


we find 


*, 
(3. 15) 8( To »t i 3 sty) ae ; 1 Y(p,t )¥(Ho» t )at ? 


tg 


5 ops, (tte 


where 


: bay Pt. 
(3-17) X(u,t) 4g i “ir | ¥(u,t)¥(u',t)at, 
0 t, 

0 


(3-18) ¥(u,t) = aie +3 ; * ai / * (us t)x(utt)e 


+ 


iplying the modified forwa: 


and e 


equation (2.20) by 
integrating over + from T) to %), we get Eqs. (3.12) — 
(3.14). Furthermore, multiply the modified forward equation 
(2.21) by _ and A ig 2 and integrate with | 
‘respect to 1 over (To»7,)- Then, we have Eqs. (3.15) — | 
(3-18). | 
In @ manner similar to that used for the forwar equations, 
‘unetions expressed in terms of | 


shall derive the S and fT. 
the and Yfunctions piper a aaa er 
Multiply (2.23) by e or"? ana ee) ° 


integrate with respect to + over (*,,1,)- 


nm we obtain 


pt. 
(3.19) S(T ost) solo) -/[ 1 X(M,t)X(Hp,t) 


(n+t)/u! | 
| d 
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(5-20) T(t 5HsHty) = a desea tty) /u oe 
A) 
where 


| (3.21) X(u,t) = 1 + f° aay : "4 X(t) x(a", 6) 


(3.22) ¥(u,t) = 


3 
X(u,t)¥(u', tle 


the other hand, mitiplying Eq. (2.27) by - og 


rating over +t from 1, to 1), | 


(5-23) S(t, tosh) = f “2 vlust)elug, tae, 
€. 
0 


(3.24) T(t), 7)s4o) = ‘ "2 (it, t)X(itpst)e © 0/0 at, 


(3-25) X(ut) = 1+ fl 


! (3.26) ¥(u,1) = Pe aA + af . We | 
6 


iT f i Y(u,t)¥(p', tat, 
* 


2 v(u,t)x(u' ere 


Vu 
dt. 


a" ‘ 


Thus, allowing for P(Ms%,7,7,) ef ‘plasty tot y—t ae) 
and the principle of reciprocity given by (3.9) and (3.10), 
we see that the laws of diffuse reflection and transmission 
based on the forward integ 
to those given by the backwaré 


odifferential equations are equal 
equations. Hence, the diffuse | 
reflection and transmission of parallel by a finite 


fy process that 
is reversible with respect to the optical depth. 

On differentiating Eqs. (3.12) and (3.15), (3.16) with 
“penpest to %, and combining them apprepriately, we obtain 
Eqs. (82), (83) in [11], Chap. VII. Eqs. (2.20) and (2.21) 
2 imilar to those provided by Sobolev [15] in the non— 


eomservative case. 
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